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NUMERICAL AND EXPERIMENTAL STUDIES OF THE NATURAL 
CONVECTION WITHIN A HORIZONTAL CYLINDER 

By Roger B. Stewart, Alexander P. Sabol, and Lillian R. Boney 
Langley Research Center 

SUMMARY 

Numerical solutions are obtained for the quasi- compressible Navier- Stokes equa- 
tions governing the time-dependent natural convection within a horizontal cylinder. The 
early flow development and wall heat transfer are obtained after a uniformly cold wall 
is imposed as a boundary condition on the cylinder. Results are also obtained for a 
time-varying cold wall as a boundary condition with windward explicit differencing used 
for the numerical solutions. The viscous truncation error associated with this scheme 

is controlled so that first-order accuracy is maintained in time and space. The results 
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encompass a range of Grashof number from 8.34 x 10 to 7 x 10 , which is within the 
laminar flow regime for gravitationally driven fluid flows. Experiments within a small- 
scale instrumented horizontal cylinder revealed the time development of the temperature 
distribution across the boundary layer as well as the decay of wall heat transfer with 
time. Agreement between temperature distributions obtained experimentally and numer- 
ically was generally good. The time decay of the dimensionless ratio of the Nusselt 
number to the one-fourth power of the Grashof number is found both numerically and 
experimentally, and good agreement is obtained between these two results over most of 
the cylinder wall. The numerical results indicate that the fluid exhibits a strong tendency 
to resist first-order motion within the inner core region. The early establishment of a 
shallow, positive, upward temperature gradient within the core enhances its stability. No 
first-order vortical motion is induced by the boundary layer and this is attributed in part 
to the fluid deceleration near the bottom of the cylinder along with expulsion of fluid 
from the boundary layer in the lower portions of the cylinder. 

INTRODUCTION 

Natural convection of fluids within closed systems has been extensively studied 
during the past 20 years. Reference 1 provides a detailed account of the work done 
prior to 1964. More recently, numerical solutions and experiments have provided 
insight into the fluid behavior for a variety of wall boundary conditions in enclosed sys- 
tems under steady-state conditions (refs. 2 to 6). 



The work concerning flow within cylindrical geometries reported in references 2, 

3, 5, and 6 has attempted to relate the dimensionless heat transfer, or Nusselt number, 
to the Grashof number, which describes the ratio of thermally driven gravitational effects 
to fluid viscous effects. A number of nonuniform wall boundary conditions have been 
investigated for the purpose of gaining insight into the coupling between the inner core 
flow and the primary gravitationally driven flow that occurs within the wall boundary 
layer. These studies have dealt exclusively with steady-state fluid flows. 

A problem closely related to the one in the present study is reported in reference 4, 
where numerical solutions were obtained for a horizontal cylinder with one-half of the 
cylinder wall heated and one-half cooled. Both the unsteady and steady -state flow were 
studied numerically, but comparisons with experiments were obtained only for the 
steady -state conditions. Very little work has been reported that compares the numerical 
and experimental development of unsteady natural convection in enclosed systems. 

The purpose of the present study was to investigate the fluid motion and wall heat 
transfer for a uniformly cold wall that was imposed as an initial condition at time zero. 
The resulting two-dimensional, unsteady laminar fluid flow is an idealization of the prob- 
lem first reported in reference 7, which considered the feasibility of developing a facility 
using stored arc-heated air to duplicate Mach 8 to 10 flight conditions. The heat loss 
to the walls of a large storage tank resulting from natural convection could conceivably 
place a size and input -power restriction on the development of such a facility. The need 
for a detailed understanding of the fluid behavior in a storage chamber motivated the 
present study. The development of the numerical solution technique is described in 
reference 8 and a stability analysis of the coupled finite-difference equations is presented. 
The present work provides numerical results describing the unsteady velocity- and 
thermal-field behavior within a horizontal cylinder and includes a quantitative assessment 
of the early fluid stratification that leads to a positive vertical temperature distribution 
in the core flow. Experimental measurements of the thermal behavior of the fluid were 
made within a small-scale horizontal cylinder. 

SYMBOLS 


A area 

c specific heat of cylinder wall material 

c specific heat at constant pressure 

r 


2 



dimensionless constant (called dimensionless heat-transfer group in ref. 4). 


w Gr 


diameter of cylinder 


gravitational acceleration 


heat-transfer coefficient 


thermal conductivity 


length 


Grashof number, 


P-w 3 T i - T v 


Nusselt number, -g- 
c d M- 

Prandtl number, -f— - 
k 

pressure 

dimensionless pressure, 
dynamic pressure 


dimensionless dynamic pressure, 


g Papj*. 


radial distance from center of cylinder 


dimensionless radial distance from center of cylinder 


dimensionless distance between two successive radial grid points 


dimensionless time, t 



dimensionless time increment 


3 




u 

u 

V 

/w 

V 
X 

a 


temperature 

temperature of outer wall surface given by equation (13) 
azimuthal velocity component 
dimensionless azimuthal velocity component, 
radial velocity component 

dimensionless radial velocity component, — 

\jZP av w 

distance from inner surface of cylinder wall toward outer wall surface 
temperature difference, T i - T w 
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A 9 
M- 


v 

P 

P 

4> 

$ 


azimuthal coordinate 

distance between two successive azimuthal grid points 
dynamic viscosity 
kinematic viscosity 
fluid density 

dimensionless fluid density, pjp^ 
dimensionless temperature function, 

dimensionless temperature function, 



Subscripts: 

i initial value at time zero 
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rain minimum 

max maximum 

w value at cylinder wall 

PHYSICAL DESCRIPTION OF THE PROBLEM 

The geometry and nomenclature of the horizontal cylinder to be studied are shown 
in figure 1, The cylinder is of infinite length to allow a region of two-dimensional flow 
to exist. When a hot gas is stored in a horizontal cylinder, appreciable azimuthal gas 
flow occurs because of the imbalance between the gravitational body forces and the 
existing hydrostatic pressure gradient in the fluid. The gas is initially at rest with a 
balance between the body forces and the hydrostatic pressure gradient and is initially at 
a uniform temperature T^. At time zero a uniformly cold wall is imposed on the cylin- 
der and the resulting conduction of heat out of the gas near the wall causes the convec- 
tive motion to begin. As the flow develops, a thin boundary layer is formed at the wall 
and thickens with time. The initially motionless inner fluid (the core) is driven by the 
boundary-layer flow and gives up energy to the heat-conducting boundary-layer fluid. 

At large values of time, the boundary layer will affect the innermost regions of the core 
flow, and eventually the gas can be expected to give up all its excess energy to the cold 
wall. Although these late values of time are not studied here, it can be anticipated that 
the fluid will be once again at rest with a uniform temperature equal to the wall 
temperature. 

Several distinct features are intuitively apparent. Because of the uniform wall 
temperature, a plane of symmetry, with the dividing line running vertically upward 
through the center of the cylinder, is immediately established. For the azimuthal veloc- 
ity component, a stagnation line coincides with the midplane of symmetry. In each half- 
cylinder there are two stagnation points which will be at the intersection of the midplane 
of symmetry and the wall. For a no-slip wall boundary condition these points will not 
appear as singularities because they are immersed in the boundary layer at its very 
limit, the wall. 

The principal motion in the boundary layer will be azimuthal and the induced motion 
will be radial. Because of the small coefficient of viscosity for air, very low velocities 
in the core may be expected. Both inward and outward radial velocities will occur across 


*A two-dimensional flow within a horizontal cylinder has been observed by Brooks 
and Ostrach (ref. 6). 
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the boundary layer. In addition, it might be anticipated that some stratification of the 
flow will occur in the lower portions of the cylinder. 

THE GOVERNING EQUATIONS 


The model chosen for this study is that of a viscous, heat-conducting, quasi- 
compressible fluid that conforms to the Boussinesq approximation. For small differ- 
ences between the gas temperature and the wall temperature, the density may be taken 
as a function of temperature only and considered as a variable only where it modifies the 
body-force terms in the equations describing conservation of momentum. This approxi- 
mation has been investigated extensively in references 3, 4, and 9 to 11. 

The quasi-compressible Navier-Stokes equations applicable to the case of large 
Grashof numbers and small differences between gas and wall temperatures, in cylin- 
drical coordinates and in nondimensional form, are 

Azimuthal momentum, 


dl 9r 


u 0u 
r 90 


s -<f> sin 9 - 


1 - 8q 

r 90 


Radial momentum, 



193 
r 9? 


u_ i_ a 2 {p 

f 2 rWy 


9r 


u 2 

- (p cos e 
r 


Mass continuity. 


Energy, 


State, 


9 (rv) . 
9f + 


9u 

90 


= 0 




d<P va <j) 2 9c/) 1 9 2 <p 

9 1 + 9r + r 90 ' Npr^Qr 9 ? 2 


( 1 ) 

( 2 ) 


(3) 


(4) 


P = 


1 

1 - j8 onp 


(5) 
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Initial conditions (at t = 0): 


{T = v = 0 " 
0 = 0 


( 6 ) 


Boundary conditions 
At r = 1.0: 

u = v — 0 

0 = 1.0 \ 

J 

On 9 = 0 and v: 

8u _ 9v _ 80 _ n 
90 “ 90 “ 90 


(7) 


( 8 ) 


The derivation of the system of equations (1) to (5) is given in reference 8 along 
with a justification for neglecting terms in both the momentum and energy equations that 
are known to be negligible for gravitationally driven fluid flows. 

Equations (1) to (5) along with the initial and boundary conditions given by 
equations (6) to (8) were solved by means of an explicit windward-differencing scheme 
that has been described in references 4, 8, 12, and 13. The formulation of the difference 
equations and a stability analysis are presented in reference 8. From the stability anal- 
ysis given in reference 8, a minimum allowable time step was determined to be 



The time step chosen for the calculations was 0.8 of the value for At given by 
equation (9). 
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The dimensionless heat transfer at the wall was obtained from 


h = 


A ( T 1 - T w) 


(where q is the heat-transfer rate) or 


, _ k OT\ 

h -'T, -TjW/r= 


w _i w 

0 = Constant 


Thus, 


(10) 


( 11 ) 


N- 


Nu 


( T i - T w) 


/9T\ 


0 = Constant 


( 12 ) 


where the derivative obtained from the difference-equation solutions was evaluated from 
the wall into the flow. Figure 2 shows a typical finite -difference grid network for one- 
half of the cylinder. Because of the symmetry about the midplane, calculations need 
only be made over half the cylinder as shown. 


NUMERICAL RESULTS 


The difference equations corresponding to the system (1) to (9) were solved for 
values of u, v, cp, and q over the entire grid after a uniformly cold wall was imposed 
as a condition at time zero. Equation (12) was solved for each time increment following 
time zero to provide dimensionless values of wall heat transfer. Three different grid 
sizes were used to study the difference approximation to the differential equations. These 
grid sizes were 31 by 51, 31 by 101, and 31 by 151, with the azimuthal grid remaining 
unchanged at a value of 31. The Control Data series 6600 computer system at the 
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Langley Research Center was used. Typical computational times were 1.02 x 10 grid 
point solutions per hour obtained for the 31 by 51 grid network. The real-time flow 
developed at a rate of about 1.25 seconds per hour of machine time for the early flow 
development. 


Input Values 

Table I lists the values of the input parameters that were used for the solutions to 
be presented. The values chosen for cases I and II resulted from the experimental condi- 
tions that were to be employed for comparison between the numerical results and measured 
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values. For cases I and II the variable Grashof numbers and wall temperatures shown 
are described later in this section. 


Constant Wall Temperature 

For the cases with constant wall temperature listed in table I, a step-function cold- 
wall condition was imposed over the entire wall at time zero. Figure 3(a) shows the 
azimuthal velocity distribution for case I and three different values of time at the azi- 
muthal location $ = 90°. The early time (f = 3.05) shows the distribution prior to the 
development of peak azimuthal velocities. By this it is meant that peak velocities at any 
location have not been achieved. The intermediate time (t - 4,57! shows the azimuthal 
velocity distribution at its peak value, and the late time (t = 14.5) shows the distributions 
after a decay in the velocity has taken place. Significant inward displacements of the 
maximum velocity are not apparent; however, the boundary layer thickens with time over 
about 20 grid point spacings for the 31 by 51 grid used. The viscous effects are trans- 
ported farther into the core of the fluid as time progresses. 

The distributions shown in figure 3 are typical of all the results obtained for a 
constant wall temperature. Figure 3(b) shows the velocity distribution near the bottom 
of the cylinder at an azimuthal angle of 23.7°. The momentum gathered by the fluid 
falling downward has both thickened the boundary layer and increased the peak azimuthal 
velocity. At values of 9 less than 23°, the fluid rapidly decelerates and comes nearly 
to rest in the lower part of the cylinder. Figure 3(c) shows the velocity distribution near 
the top of the cylinder. The boundary layer is well defined, but the azimuthal velocities 
are substantially lower than in the bottom portion of the cylinder near the wall. 

Figure 4(a) shows the radial velocity distribution for case I of table I. Here it is 
seen that the peak velocities are an order of magnitude less than corresponding azi- 
muthal velocities. The radial velocities peak at the later time (f = 14.5) in contrast to 
the azimuthal velocities, which peak near t =4.57. The negative radial velocity in 
figure 4(b) at 9 = 23.7° indicates an expulsion of fluid from the boundary layer. As 
the lower portion of the cylinder "fills" with fluid that has moved downward in the bound- 
ary layer, negative or inward radial flow takes place to accommodate the added fluid. 
Thus fluid is slowly forced into the boundary layer near the top of the cylinder because of 
displacement effects from below. Figure 4(c) shows the radial velocity distribution in 
the upper portion of the cylinder. At 9 = 161° the flow is still developing, with peak 
velocities occurring near t = 14.5. 

Figure 5 shows two radial velocity distributions within the core flow as a function 
of azimuthal location. The small values of v over the entire azimuthal grid for both 
early and late times are representative of the core resistance to induced fluid motion. 
The azimuthal velocity in the boundary layer is many orders of magnitude larger than the 
radial velocity in the core. 
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The general picture of the flow field involves three major features: (1) a boundary- 
layer development near the wall due primarily to downward azimuthal fluid flow near the 
wall, (2) an induced radial velocity that feeds fluid into the developing boundary layer in 
the upper and middle azimuthal locations of the cylinder and ejects fluid out of the bound- 
ary layer at lower azimuthal locations, and (3) a core region that strongly resists first- 
order motions and only very slowly forces fluid at lower levels to rise and enter the 
boundary layer. This is a striking example of a flow in which the principal motion is 
confined to the boundary layer. 


The dimensionless temperature <f> is shown as a function of radial distance in 
figure 6 for three different azimuthal locations. Here the thermal boundary layer is less 
well developed than the velocity layer at a given time. The conduction terms in the 
energy equation dominate the early development of the thermal boundary layer because of 
the very low fluid velocities that occur shortly after the cold-wall initial condition is 
imposed. 


The early behavior of wall heat transfer is of principal interest, and figure 7 shows 
the Nusselt number decay with time for four different azimuthal locations. These 
results appear to be a valid physical representation of the flow development because the 
convective terms in the energy equation were found to be negligible when very low fluid 
velocities are present in the early development. The result is an energy balance that is 
controlled by conduction out of the fluid to the wall until the velocity field develops. The 
increasing value of Nusselt number with increasing 9 gives a clear picture of the posi- 
tive upward temperature gradient within the boundary layer. This positive upward gra- 
dient also exists in the core fluid, as indicated by figure 8. The early establishment of 
this upward gradient produces a thermally stable core which tends to resist downward 
motion. The low velocities calculated for the core flow are, in part, a result of this 
thermal behavior. Figure 9 shows a plot of N /Ncr^ 4 as a function of time - The 
cold fluid entering and residing in the lower portion of the cylinder causes a rapid and 
nearly linear decay of the dimensionless heat-transfer function. For 9 near 90° and 
up to 180° the function decays less rapidly because of the flow of warm fluid into the 
boundary layer from the core. At very large times the heat-transfer function must 
asymptotically approach zero because the fluid gives up all its excess energy to the wall. 


It is apparent that c w = 


~nr -® as t - 0 

N 1/4 
w Gr 


imposed as an initial condition. 


because of the step-function cold wall 


The calculated value of at early times near t = 0 is dependent on the radial 
grid spacing used, as might be expected. A comparison of the effects of grid refinement 
is shown in figure 10. In figure 10(a) the value of T/t^ as a function of radial distance 
is shown for three different grid networks, where the radial coordinate grid spacings are 


10 



0.02, 0.01, and 0.0067. The azimuthal location Q is 90° and the time t is 2.7. As 
might be expected, a finer grid network steepens the temperature gradient near the wall. 
The convergence toward a solution that is independent of grid size is indicated by fig- 
ure 10. Of interest is the velocity distribution shown in figure 10(b). Here, also, con- 
vergence is indicated by the results for the refined grid. The velocity peak for the 
coarse grid, however, is above the peaks for the finer grids. The velocity distribution 
near the wall appears to be adjusting itself from an "overshoot'’ where Ar = 0,02 to a 
convergent result as Ar becomes smaller. Similar results have been computed for 
Grashof numbers of 8 x 10 4 and 7 x 10 7 for the grid spacings given above. 

Case IH of table I was computed for a Grashof number of 8.34 x 10 4 . The computa- 
tional time increases considerably at this relatively low Grashof number. Figure 11(a) 
shows representative azimuthal velocity distributions for this case. These results for 
the lower Grashof number can be interpreted as being due to a more viscous fluid flow. 

At a dimensionless time t of 3.07 the boundary layer is considerably thicker than the 
comparable boundary layer of case I (fig. 3(a)). Even at the early time of t = 1.31 the 
viscous effects are more pronounced than for the higher Grashof number. Figure 11(b) 
shows typical radial velocity distributions for the low Grashof number (case in). At 
t = 3.07 the distribution shows the effects of higher viscosity in that the profile is 
thicker than for the case I solution. Finally, figure 11(c) shows the distribution of the 
dimensionless temperature function. Again, the more viscous fluid of case III shows a 
thicker thermal boundary layer than that of case I. All the profiles for case IH have a 
qualitative resemblance to those for case I. The solutions for the Nusselt-Grashof 
relation at the lower Grashof number are shown in figure 11(d). These curves should be 
compared with the results shown in figure 9. 

The computations for the dimensionless dynamic pressure q indicate that the 
dynamic pressure gradient makes only a negligible contribution to the momentum balance 
within the fluid; thus, calculations for q were discontinued. This occurrence is consis- 
tent with the results reported by Ostrach in reference 1. 
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The results of the computations for case II at a Grashof number of 1.3 x 10 are 
shown in figure 12. The velocity and temperature profiles represent a result that is 
intermediate to the solutions for case IH and case I. 

Time-Varying Wall Temperature 

In order to compare the numerical solution with the experimental results, it was 
necessary to use a time -dependent wall temperature as a boundary condition. The experi 
mental wall-temperature decay on the inner face can be described by a solution of the one 
dimensional heat-conduction equation. For a wall of thickness x Q initially at a tempera 
ture 4>(x,0) = 1, where 
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( 13 ) 






the temperature history at the inner face x = 0 is given by 

,1A 


4.(0, t)= 1 - £ (-l n ) 


n=0 


erfc (pc)V2 

2yfcx 0 t \ w/ 


(14) 


With the solution to equation (14) the time dependence of the coefficient that 

appears in equations (1) and (4) may be calculated: 


1 1 



This value is calculated at each time step, and the difference equations are solved by 
using each new value of Grashof number which results from the time-varying wall 
temperature. As an example of the solution given by equation (14), figure 13 shows the 
time-dependent dimensionless wall-temperature decay. Several different limits on the 
number of terms taken in the summation of equation (14) showed rapid convergence for 
values of n greater than about 6. 

A fundamental difficulty is encountered when the error function solution given by 
equation (14) is used. At very early times the function is so close to unity that a near 
singularity is introduced into equation (15). The difference scheme is limited to small 
time steps, however, for solution of the momentum and energy equations. The wall- 
temperature decay obtained from equation (14) is numerically incompatible with the 
difference scheme and does not provide a usable wall boundary condition. This difficulty 
can be overcome by approximating equation (14) at early times by linear functions of the 
form 


$(o,t) = 0.99 - 0.050? (16) 

Figure 14 shows the azimuthal velocity distribution for case I when the wall temperature 
decays according to equation (16) rearranged: 



The solutions for a time-dependent wall temperature exhibit several interesting 
differences from those for a constant wall temperature. In effect, a boundary condition 
of a slowly cooled wall produces an early driving force that is quite small. This is 
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characterized by the time development of the Grashof number. Figure 1^6 shows a typical 
time history of the Grashof number for the input values of case I and a wall that follows a 
temperature history given by equations (16) and (17). The development of the flow field 
is directly related to the time -dependent Grashof number. The immediate consequence of 
a slowly cooled wall should appear as a flow that is less fully developed at any given time 
than the flow for which a step-function change in wall temperature has been imposed. The 
early behavior should be closely related to a flow field having a lower Grashof number. 

The boundary-layer development follows the rising Grashof number with a behavior that 
is different from that for a constant wall temperature. In effect, the driving force 
increases with time, and thus the flow experiences a longer and more pronounced 
acceleration. 

Figure 16(a) shows typical temperature distributions for the boundary conditions 
of a time -dependent wall temperature and a constant wall temperature. The details 
revealed by these curves are physically valid in that the distribution for a constant wall 
temperature, of necessity, must have a steeper slope near the wall in addition to a 
thinner thermal boundary layer due to the higher Grashof number at equal values of time. 
These features are clearly evident in figure 16(a). 

A comparison of the distributions for constant and variable wall temperatures at 
equal values of the Grashof number reveals the fundamental differences between the two 
cases. Figure 16(b) shows that even for equal values of Grashof number the variable 
wall temperature still has the effect of producing a more viscous flow throughout the 
entire time development than the constant wall temperature. A comparison with the 
results for a constant Grashof number is difficult because of the fundamental difference 
between the two case histories. When the two boundary-layer thicknesses are compared 
at equal values of r and time, the results shown in figure 17(a) are obtained. The 
difference in the peak values of the dimensional velocity is of course even greater than 
that shown in this figure because of the different value of reference velocity used in the 
nondimensionalization, that is, u = \Jgj3ar^ u. A comparison of the dimensionless 
azimuthal velocity distributions at equal values of both Grashof number and time is shown 
in figure 17(b). As expected, the distribution for the time-varying wall temperature 
shows a thicker profile but a lower peak velocity. Similar results were observed for the 
temperature distribution, which is also an indirect result of the viscous behavior of the 
fluid. 

Perhaps the most significant single consideration regarding the energy transfer to 

the walls is the time -dependent Nusselt-Grashof correlation. Heliums (ref. 4) estab- 

1/4 

lished a constant in the relation N Nu = c w N Q r • For one particular steady flow within 
a horizontal cylinder, c w = 0.326. If a fully developed flow occurs in the present problem, 
the time -dependent values of c (T) should approach the steady value prior to a decay in 
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the fluid motion. As time increases indefinitely in the present unsteady flow, c w (T) -0. 
The fluid will then have given up all its excess energy to the cold walls. Figure 18 shows 
the time -dependent behavior for the Nusselt-Grashof correlation. The relation 
N Nu/ N Gr 4 is seen t0 a PP roach the steady-state value within 1 second (t = 2.8) after 
flow initiation. The unsteady-flow value of this relation for the present problem must be 
asymptotic with zero at large values of time. This behavior is implied by the numerical 
results represented in figure 18. 

Figure 19 illustrates the time history of "fluid particles" within the flow field by a 
sequence of photographs that were taken from an oscilloscope display of the particle dis- 
placements computed from the numerical solutions. As shown in these figures, particles 
were located along rays separated by an azimuthal distance of w/8 at t = 0. From 
r = 0 to r = 0.64 the radial spacing between particles Ar was 0.04. From r = 0.64 
to r = 1.0 the radial spacing Ar was 0.01. A two-dimensional interpolation was used 
to obtain particle displacements between grid locations for the 31 by 101 mesh in these 
calculations. The time-dependent wall -temperature decay for case I of table I was used. 
The sequence of photographs shows the fluid deceleration and stratification near the 
bottom of the cylinder. Regions in which the number of particles increase represent 
higher fluid density, as might be expected. Figure 19(c) is a plot of the displacements 
at 6.6 seconds. In this figure the particles belonging to each original ray have been 
faired in to illustrate the displacement profiles. Real-time motion pictures were made 
from plots such as these to give a physical picture of the flow-field development. Because 
of the low velocities within the core flow, large values of time must be obtained with the 
numerical solutions before substantial particle displacements can be observed within the 
core. 


EXPERIMENTAL APPARATUS AND PROCEDURE 

The instrumented stainless-steel cylinder used for measuring temperature distri- 
butions in the gravitationally driven flow field is shown in figure 20. Experiments were 
made with the cylinder in a horizontal position, with rotation of the cylinder allowed in 
order to position the thermocouple probes at different locations. The inside diameter of 
the cylinder was 15.1 cm and the wall thickness was 0.64 cm. The cylinder was 152 cm 
in length with stainless-steel caps welded on each end. An optical transit was used to 
provide precise orientation of the cylinder with respect to the Earth's gravitational field 
vector. Eighteen thermocouple probes were installed over the entire length of the cylin- 
der to provide checks on any end-wall effects that might violate the assumption of two- 
dimensional flow. Measurable end-wall effects, signifying three-dimensional flow, 
extended less than 2 cylinder diameters from each end wall. The cylinder was also 
instrumented with a series of 30-gage copper- constantan thermocouple probes spaced at 
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i/d = 5 to allow temperature measurements to be made across the boundary layer and 
into the core. The probes were installed as shown in figure 21 with the intent of intro- 
ducing a minimum of disturbance to the flow passing over the probe sensor regions. At 
large times the probes could affect the flow recirculating from the core into the boundary 
layer, but the very low fluid velocities and core stratification found in the numerical solu- 
tions indicated that such disturbances should be negligible. The errors associated with 
conduction, radiation, thermal lag, and velocity losses were investigated in reference 8. 
There it was found that for the test conditions reported, the radiation and thermal-lag 
errors are significant, whereas conduction and velocity errors remain negligible. 

The outer jacket surrounding the cylinder, the hoppers, and the manifold system are 
shown schematically in figure 22. The 2.54-cm-thick Micarta rings at the bottom of the 
hoppers served to insulate the outer jacket from the cold liquid being stored in the hop- 
pers. At time zero the plug valves were opened manually and the liquid in the hoppers 
was dumped around the cylinder walls. The dumping typically took less than 1 second. 

It was not possible to provide a step-function input to the inner cylinder wall because of 
the finite thickness of the wall. As mentioned in the previous section, the one -dimensional 
heat- conduction equation was applicable to a description of the decay of the inner-wall 
temperature. Once the cold liquid was discharged from the hoppers onto the cylinder 
wall, a circulation pump was activated to destroy temperature gradients in the cold 
liquid. In particular, a thermal boundary layer is established between the cylinder wall 
and the liquid, and this would preclude an accurate prediction of the inner-wall tempera- 
ture history. 

Considerable care was taken in locating the boundary -layer thermocouple probes. 
These fine wires were quite easily displaced from a fixed position, so that after locating 
each wire in the correct position, no further instruments or probes were inserted inside 
the cylinder. Also, charging and purging of dry air in the cylinder were done at a slow 
rate, so that convection velocities remained below 61 cm/ sec at all times. After the test 
gas was introduced into the cylinder, a waiting period of 2 hours followed to insure that 
all motion had damped out within the cylinder. From the junction box the signals were 
transferred to an analog -to- digital data recording system and then read out on computer 
tape. 

One method for varying the Grashof number from one test to another was accom- 
plished by changing the pressure within the cylinder. Each test was repeated a minimum 
of three times, and calibrated pressure transducers were used to monitor the overall 
tank pressure during each test. The j transducers had a ± 5-percent error limitation, but 
for all duplicate tests they showed a repeatability better than 5 percent. The coolant 
temperature was monitored by a thermocouple located in the annular space between the 
manifold and cylinder wall. The temperature of the outer cylinder wall was also 
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monitored with a thermocouple installed on the surface. The experimental data reported 
here were obtained under conditions that produced a minimum ratio of wall temperature 


to gas temperature 


(KL ■ M36 ) 


when a nearly steady wail temperature had been 


achieved after the application of the initial condition. 


COMPARISON OF NUMERICAL AND EXPERIMENTAL RESULTS 


Comparison of the numerical and experimental results necessitated a consistent 
formulation of the dimensionless variables as well as identical wall boundary conditions. 
Also, the numerical solutions must have the initial conditions imposed at the same time 
zero that actually occurred in the experiments. It was found that within 0.4 second 
after application of the cold liquid to the cylinder wall, temperature drops were detected 
within the gas near the wall. Thus, time zero could be determined within 0.4 second by 
thermocouple signals alone. For comparison with the numerical results the heat- 
conduction equation was used to determine time zero for the experiments by allowing 
time zero for the conduction problem to occur when the liquid was first dumped from the 
hoppers. Time zero for the numerical computations and the fluid was taken when a 0.11K 
decrease in wall temperature had occurred as a result of conduction. 

The foregoing considerations were used to obtain the measured and computed 
temperature distributions for an azimuthal angle of 90° shown in figure 23. These results 
correspond to case I of table I with the exception of a variable wall temperature, as men- 
tioned above. The numerical solutions for 1.2, 4.0, and 8.0 seconds of real time com- 
pare favorably with the measured results. The large computing time required for the 
program prohibited carrying out the numerical solutions to larger values of time. The 
gradual thickening of the thermal boundary layer is seen in figures 23(a) to 23(c). 

Because the core acts as a reservoir of warm fluid for supplying the boundary layer, the 
temperature distributions within the layer tend to retain their profiles over a long period 
of time after the velocity field starts decaying. 

To illustrate the temperature decay, figure 24 shows typical plots of both the 
experimental and theoretical temperature at selected azimuthal locations and at the 
closest radial location to the wall. The thermal behavior shown in figure 24(a) clearly 
indicates the inflection produced by the flow as the velocity distribution shifts toward a 
more developed profile. Similar behavior is seen in figure 24(b) at 9 = 70°. Finally, 
at the bottom of the cylinder, where the azimuthal motion of the fluid ceases, the tempera- 
ture decay shown in figure 24(c) appears. 

The numerical results predict a more rapid decay than was actually measured. 
Several reasons for this discrepancy are possible. The most likely source of error lies 
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within the framework of the differencing scheme near the midplane of symmetry. Because 
both the azimuthal and radial velocities are negative near the bottom of the cylinder, for- 
ward differencing is used throughout the momentum and energy equations. But forward 
differencing carries the finite-difference calculations into regions of largest change and 
tends not to balance with the backward grid points that for the present physical system lie 
in regions of lesser change. Near the midplane the azimuthal velocity goes to zero, but 
forward differencing tends to override this, partly because of the finite grid size; thus it 
appears that errors may be largest near the bottom midplane of symmetry, In this 
respect the solutions are dependent upon the cylindrical geometry being considered. 

This dependence represents a limitation, or at least an undesirable aspect, of one-sided 
finite -difference techniques. 

It is of some interest to observe the experimentally measured pressure decay within 
the cylinder. This decay represents a three-dimensional (closed-volume) phenomenon 
which is in a sense primarily dependent on the two-dimensional flow field that transports 
heat to the cylinder walls. Figure 25 shows a typical measured and calculated pressure 
decay for a maximum Grashof number of 7 x 10 7 (case I). The pressure decay is a 
direct measure of the overall energy loss from the fluid. The average density in the 
enclosed cylinder remains constant, and thus values of the average tank pressure are 
directly proportional to the average fluid temperature within the cylinder. 

The single most useful comparison made between the numerical results and the 
experiments is that of the Nusselt-Grashof relation for the time-dependent flow. Fig- 
ure 26 shows the measured time-dependent variation of this dimensionless quantity. The 
comparison is favorable. The Nusselt-Grashof relation approaches the steady-state value 
within 4 seconds (t = 10) after flow initiation and is asymptotic with zero at large values 
of time. 

If figure 26 is compared with figure 18, the effect of a steadily decaying wall tem- 
perature is apparent. The curve shown in figure 26 results because lower values of heat 
transfer occur when the wall is cooled unsteadily than when a step-function cold wall is 
applied, such as in the conditions for figure 18. 

Examples of typical experimental radial temperature distributions are shown in 
figure 27. In figure 27(a) at t = 4.0 seconds, the wall has cooled to about 0.983 of the 
initial temperature, and thus the distribution terminates at the wall (r = 1.0J in the 
manner shown. 

The thermocouple located at r = 0.2 was considered essential for monitoring the 
core temperature. If a very low ratio of T^Tj had been utilized, it could be anticipated 
that fluid would fall downward from the uppermost walls of the cylinder. In such a case 
a negative vertical temperature distribution within the core might have been observed. 

For the value of T w /Tj of these experiments the flow was almost completely confined 
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within the wall boundary layers, and no negative upward distributions were measured. 
Figure 27(b) shows the temperature distribution at 8.0 seconds after flow initiation. The 
measured profile has broadened more than the numerical solutions predict. The agree- 
ment is reasonable, however, considering the high degree of nonlinearity of the problem 
when the walls are cooled unsteadily. Figures 28(a) and 28(b) show a similar result for 
6 = 70°. 

5 

A series of experiments at Grashof numbers up to about 6.5 x 10 were made, and 
some results are shown in figure 29. The very slow time development of the flow field 
is indicated both by the numerical calculations, which take a computational time almost 
an order of magnitude longer than the case I solution, and by the experiments. The agree- 
ment between the numerical solutions and the experiments is not as good as that for the 
results at the higher Grashof number. The very slow development of the flow reduces the 
accuracy at early times because of the small differences in temperature that must be 
measured. Figure 30 is a plot of the time history of the boundary-layer temperature at 
r = 0.966. The slow decay is an example of the behavior at the lower Grashof number. 
However, the more rapid reduction of the ratio of gas temperature to wall temperature 
for the lower Grashof number produces an undershoot of the Nusselt-Grashof relation. 
This is due to the differences between the velocity-field development and the thermal- 
field development and is shown in figure 31. An initial response to the near singularity 
at early times due to low values of the Grashof number causes the undershoot in the 
numerical results. Thus, although the disturbance is damped by the difference scheme, 
the accuracy of the results is severely lessened by the destabilizing influence of a very 
low Grashof number at early times. Restated, the computations for the low Grashof 
number are characterized by a relatively rapid wall cooling which raises the wall heat 
transfer more rapidly than the one-fourth power of the Grashof number. The Nusselt- 
Grashof behavior in both figures 26 and 31 must be viewed as being representative of one 
specific flow configuration and also a very specific boundary condition. The difference 
between the results shown in figures 26 and 31 indicates that a single correlation equation 
from the numerical results presented here may not be possible to achieve for the entire 
range of unsteady laminar natural convection within a horizontal cylinder. 

CONCLUDING REMARKS 

The present investigation has revealed the dynamic and thermal behavior of a con- 
fined fluid subjected to gravitational body forces. For an infinite horizontal cylinder with 
uniformly cold walls imposed as an initial condition, the principal features of the flow 
involve a rapid development of the boundary layer adjacent to the walls, a tendency for 
fluid stratification in the lower regions of the cylinder, and a slow decay of the velocity 
and thermal fields with time. Experiments made in this investigation substantiate the 
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thermal behavior predicted by numerical solutions to the quasi-compressible Navier- 
Stokes equations for a time -dependent wall-temperature decay. Because the thermal and 
velocity fields are strongly coupled, the experimental findings imply that the velocity field 
may also be accurately described by the numerical results. The following remarks can 
be made from this study: 

The relationship between the Nusselt and Grashof numbers is found both numerically 
and experimentally for unsteady natural convection within a horizontal cylinder subjected 
to uniformly cold walls as boundary conditions. For the condition of a constant cold wall 
from time zero, the numerical results show that the Nusselt-Grashof relation approaches 
the steady- state value about 1 second after commencement of the flow down the cylin- 
der wall. For the boundary condition of a time-dependent wall temperature, the 
Nusselt-Grashof relation approaches the steady-state value about 4 seconds after the 
commencement of flow down the cylinder walls. 

The combination of early establishment of a positive upward temperature gradient 
within the central core flow and the effects of fluid viscosity produces a strong resis- 
tance to induced fluid motion within the core. After initiation, the flow within the bound- 
ary layer develops rapidly, and decays over a long period of time. 

The dynamic pressure gradient terms in the Navier-Stokes equations for azimuthal 
and radial momentum are negligible for the conditions considered in this investigation. 

For natural convection with very low velocity, windward finite differencing, which 
gives first-order accuracy in time and space, is a suitable numerical scheme when posi- 
tive and negative velocities occur. The large computing time required for such flow 
fields rules out more complex differencing schemes at the present time. 

The differencing scheme appeared to lose accuracy for the flow near the bottom of 
the cylinder, where both the radial and azimuthal velocities are negative and the azimuthal 
velocity decelerates to zero at an azimuthal angle of zero. 

No first-order motion was found within the core flow. The absence of this motion 
is attributed to the boundary conditions, which require azimuthal deceleration of the 
boundary -layer flow near the bottom of the cylinder. Both the theoretical and experimen- 
tal models establish a midplane of symmetry that satisfies the boundary conditions imposed 
in this study. 

The present studies are applicable strictly to the regime of laminar natural con- 
vection. The limit of laminar flow over cylinders in the external flow situation has been 
observed to occur at Grashof numbers near 3 x 10 9 . For the high-pressure full-scale 
facility discussed in AIAA Paper No. 71-254, the Grashof number near the end of the 
charging cycle would be slightly in excess of this value, and turbulent heat-transfer rates 
could be expected. Experimental steady-state heat-transfer results, however, show only 
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a slight difference in slope due to the change from laminar to turbulent flow. In addition, 
the time history of the rate of increase in heat transfer in the full-scale charging cycle 
would produce a time-dependent Grashof number even for a steady-state wall temperature. 

Finally, a full-scale facility could introduce axial flow within the chamber as a 
result of the method of injection of gas during charging. Only if the axial velocities were 
much less than the azimuthal velocities induced by a cold-wall boundary condition, would 
the overall heat transfer be expected to show small differences from the value for a 
truly two-dimensional flow field. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 24, 1974. 
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TABLE I.- PARAMETER VARIATIONS FOR NUMERICAL STUDIES 


Case I: 

Grashof number 

Prandtl number 

T w/ T i 

Tw 

r w , cm 

P t , atm* 

t, sec 

Case II: 

Grashof number 

Prandtl number 

TwAi 

Tw 

r w , cm 

atm 

t, sec 

Case HI: 

Grashof number 

Prandtl number 

Tw/Ti 

T w 

r w , cm 

Pj, atm 

t, sec 

*1 atm = 101.3 kN/m2. 


7 x 10? 

0 to 7 X 10? 

0.715 

0.715 

0.936 

Variable 

Constant 

*(°> T X T 1 - T xo )+ T xo 

7.6 

7.6 

6.74 

6.74 

0.35 1 

Variable t 

1.36 X 10 6 

0 to 1.36 x 106 

0.715 

0.715 

0.936 

Variable 

Constant 

*(0,T)(Ti-t X o) + t xo 

7.6 

7.6 

1 

1 

0.35t 

Variable t 


8.34 x 104 
0.715 
0.936 
Constant 
3.1 
1 

0.22 1 
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Figure 1.- Geometry of horizontal cylinder. 










(a) 9 = 90°. 

Figure 4.- Radial velocity variation with radial distance. Case I; 
N Gr = 7 x 10 7 ; grid, 31 by 51. 























Radial distance, r 

(a) Temperature distribution. 


Radial distance, r 

(b) Velocity distribution. 


Figure 10.- Effects of radial grid refinement for a Grashof number of 1,36 x 10°. 

Case Hr t = 2.7. 
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(a) Azimuthal velocity distribution. Case HI; NQ r = 8.34 x 10^ 

9 = 90°; grid, 31 by 51. 

Figure 11.- Results for natural convection at low Grashof number 
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(b) Radial velocity distribution. Case III; NQ r = 8.34 x 10^; g = 90°; grid, 31 by 51. 
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Radial distance, r 

(c) Temperature distribution. Case III; Nqj, = 8.34 x 10^; 9 = 90°; grid, 31 by 51 

Figure 11.- Continued. 















Dimensionless time, t 


Figure 15.- Time development of the Grashof number for unsteady wall temperature 

Case I; grid, 31 by 51. 
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(a) Variable Grashof number. 

Figure 16.- Temperature distributions for a constant and a time -dependent wall 

temperature. Case I; 6 = 90°. 
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(b) Equal Grashof numbers = 

Figure 16.- Concluded. 













t = 0; N Gr = 0 



(a) t = 0 and 3.86 sec. 

Figure 19.- Oscilloscope photographs of fluid-element displacements from the numerical 
solutions for natural convection within a horizontal cylinder. Time-dependent wall 
temperature; case I. 


t = 5.13 sec; N Gr = 3.6 x 10 



t= 6.15 sec; N Gf =4.9xl0 7 

(b) t = 5.13 and 6.15 sec. 
Figure 19.- Continued. 
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Figure 21.- Locations of boundary -layer and core thermocouples within horizontal cylinder. 
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Figure 22.- Schematic diagram of apparatus for natural convection experiments. 




Radial distance, r 
(a) t = 1.2 sec. 

Figure 23.- Temperature distribution for a time-dependent wall temperature 

Case I; e = 90°. 
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Figure 25.- Measured pressure decay due to natural convection within a horizontal 

cylinder. Case I; (hU \ = 7 x 10 7 . 
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Figure 26.- Nusselt-Grashof relation for natural convection within a horizontal cylinder. 

Case I; (N G ) = 7xl0 7 ; 9 = 90°. 
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(a) t = 4.0 sec. (b) t = 8.0 sec. 

Figure 28.- Radial temperature distribution at 9 - 70°. Case I; = 6 x 10?. 
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Figure 30.- Decay of boundary -layer temperature at 6 = 90°, r = 0.966, and 

(Nrjr'l = 1.3 x 10 6 . Case II. 
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Figure 31.- Nusselt-Grashof relation for 0=90° and (NrJ = 6.5 x 10 5 . Case n. 
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